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Abstract. We present a complete analytic solution to the problem of count- 
ing planar graphs. We prove an estimate g n ~ g-n~ 7 / 2 7 n n! for the number g n 
of labelled planar graphs on n vertices, where 7 and g are explicit computable 
constants. We show that the number of edges in random planar graphs is 
asymptotically normal with linear mean and variance and, as a consequence, 
the number of edges is sharply concentrated around its expected value. More- 
over we prove an estimate g(q) ■ n~ 4 ' j y(q) n n\ for the number of planar graphs 
with n vertices and [qn\ edges, where -f(q) is an analytic function of q. We also 
show that the number of connected components in a random planar graph is 
distributed asymptotically as a shifted Poisson law 1 + P(u), where v is an ex- 
plicit constant. Additional Gaussian and Poisson limit laws for random planar 
graphs are derived. The proofs are based on singularity analysis of generating 
functions and on perturbation of singularities. 



1. Introduction and statement of results 

In this paper we obtain a precise asymptotic estimate for the number of labelled 
planar graphs on n vertices, and we establish limit laws for several parameters in 
random labelled planar graphs. In particular, we show that the number of edges in 
random planar graphs is asymptotically normal, and that the number of connected 
components in a random planar graph is distributed asymptotically as a shifted 
Poisson law. Additional Gaussian and Poisson limit laws for random planar graphs 
are derived. 

From now on, unless stated otherwise, all graphs are labelled. Recall that a 
graph is planar if it admits an embedding in the sphere. We remark that we 
consider planar graphs as combinatorial objects, without referring to a particular 
topological embedding. 

Let g n be the number of planar graphs on n vertices. A superadditivity argu- 
ment ^21 shows that the following limit exists: 

7 = lim (g n /n\) 1/n . 

n — >oo 

Until recently, the constant 7 was known only within certain bounds, namely 

26.18 < 7 < 30.06. 



Date: February 1, 2008. 

1991 Mathematics Subject Classification. Primary 05A16, 05C30; Secondary 05C10. 
Key words and phrases. Planar graph, random graph, asymptotic enumeration, limit law, 
normal law, analytic combinatorics. 

Research supported in part by Projects BFM2001-2340 and MTM2004-01728. 

1 



2 



OMER GIMENEZ AND MARC NOY 



The lower bound results from the work of Bender, Gao and Wormald jT]. They 
show that, if b n is the number of 2-connected planar graphs, then 

lim (b n /n\) 1/n » 26.18. 

n — >oo 

Hence 7 is at least this value. 

The upper bound is based on the fact that an unlabelled planar graph on n 
vertices can be encoded with at most cm bits for some constant a. If this is the case 
then g n < 2 an n\, and so 7 < 2". The first such result was obtained by Turan [T?)] 
with the value a — 12. This has been improved over the years and presently the 
best result is a w 4.91, obtained by Bonichon et al. Since 2 491 30.06, the 
upper bound follows. 

Recently the present authors ^U] were able to obtain, using numerical methods, 
the approximation 7 rj 27.2268. In this paper we determine 7 exactly as an analytic 
expression. Moreover, we find a precise asymptotic estimate for the number of 
planar graphs. 

Theorem 1. Let g n be the number of planar graphs on n vertices. Then 

(1.1) g n ~ g ■ n- 7 /V»!, 

where g w 0.4260938569 • 10~ 5 and 7 « 27.2268777685. 

The constants in the last statement are completely determined as analytic ex- 
pressions in terms of elementary functions. The proof of Thoerem^ together with 
the expressions given in the appendix, contain all the necessary details for deter- 
minining the constants. This also applies to all the remaining constants that appear 
in the paper. 

As we show later, for the number c n of connected planar graphs on n vertices, 
we have the estimate 

c n ~ c • n~ 7/2 j n n\, 
where 7 is as before and c w 0.4104361100 • 10~ 5 . 

The proof of Theorem^is based on singularity analysis of generating functions; 
see EJ - Let g n , c n and b n be as before. As we show in the next section, there are 
two equations linking the exponential generating functions 

B(x) = b n x n /nl, C(x) = CnX n /n\, G(x) = £ g n x n /nl 

The dominant singularity of B(x) was determined in p^; we are able to obtain the 
dominant singularities of C(x) and G(x), which are both equal to p — 7 _1 . 

In Section we review the preliminaries needed for the proof. In Section 
we find an explicit expression for the generating function B(x,y) of 2-connected 
planar graphs counted according to the number of vertices and edges. This is a 
key technical result in the paper, which allows us to obtain a full bivariate singular 
expansion of B{x, y). The explicit expression obtained for the function (3{x, y, z, w) 
in the statement of Lemma El suggests that we are in fact integrating a rational 
function. This is indeed the case as we explain later. 

In Section 0] we determine expansions of C(x) and G(x) of square-root type at 
the dominant singularity p, and then we apply "transfer theorems" [SjE] to obtain 
estimates for c„ and g n . 

The singular expansions of C(x) and G{x) can be extended to the corresponding 
bivariate generating functions C(x,y) and G(x,y) near y = 1. This allows us to 
prove in Section using perturbation of singularities 6 , a normal limit law for 
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the number of edges in random planar graphs. To our knowledge, this problem was 
first posed in ^Q. 

Throughout this paper, we say that a sequence of random variables X n with mean 
[i n and variance a\ has a normal limit law if the normalized variables X* = (X n — 
A*n)/°Vi converge in law to the standard normal distribution 7V(0, 1); convergence 
in law means, as usual, point-wise convergence of the corresponding distribution 
functions. 

Theorem 2. Let X n denote the number of edges in a random planar graph with 
n vertices. Then X n is asymptotically normal and the mean [i n and variance <j„ 
satisfy 

(1.2) (in ~ ku, ofj ~ Xn, 

where k w 2.2132652385 and X w 0.4303471697. The same is true, with the same 
constants, for connected random planar graphs. 

As a consequence, since a n = o(/i n ), the number of edges is concentrated around 
its expected value; that is, for every e > we have 

Prob{|X n — nn\ > en} — > 0, as n — > oo. 

Previously it had been proved that Prob{X„ < an} — » and Prob{X„ > (3n} — > 0, 
as n — > oo, for some constants a and (3. The best values achieved so far were 
a w 1.85 (shown in improving upon 0]) and /3 w 2.44 (shown in [3], improving 
upon ^3] ). Theorem |3 shows that in fact there is only one constant that matters, 
namely k. 

The previous theorem shows convergence in distribution to the normal law. How- 
ever, in this setting it is often the case that one can also prove a local limit law, 
that is convergence to the density function of the normal law. We prove such a 
local limit law and we derive large deviation estimates for the number of edges in 
random planar graphs. In the next statement, as later in the paper, g n ^ q and c n>q 
denote respectively the number of planar graphs with n vertices and q edges; p(y), 
G§{y) and C§{y) are computable analytic functions to be introduced later. 

Theorem 3. Let /i be a fixed ratio in the open interval (1,3). Take u > such 
that —up'(u)/p(u) = \i. Then, as n goes to oo, 

(1-3) fln.LjunJ ~ n\G 5 (u)- 



/2~^r(-5/2)cro 7 /2' 
where 



lP "{u) p'(u) , 2 p'{uf 



2 2/ 

a = —u — — u — -— + u 



p{u) " p(u) ' " p(u) 2 ' 
The same is true for the number of connected planar graphs c n ; \u n \ if we replace 
G 5 (u) by C 5 {u) in fO|) . 

The previous result makes more precise a recent result from 9 , where the authors 
show that 

lim ilog^lpi =a(m)! 
n— >oo n n\ 

where A(/i) is a continuous function of /i. A direct consequence of Theorem [21 is 
that 

A(/i) = -/ilog(it) - log(p(u)) 
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Figure 1 . The growth ratio of planar graphs with n vertices and 
\jin\ edges 

where u depends on fi as in the statement of the theorem. Notice that A(/i) is an 
analytic function of fi. Figure 2] shows the plot of exp(A(/z)), that is, the growth 
ratio of planar graphs with n vertices and [^"J edges. The limit of exp(A(/z)) as 
fi — ► 1 is equal to e, which is the growth ratio of labelled trees; the limit as /i — ► 3 
is equal to 256/27, which is the growth ratio of triangulations [T7| . (Tutte's result 
is for unlabcllcd triangulations, but a triangulation has at most a linear number of 
automorphisms . ) 

Next we turn our attention to the following problem, considered in 12]. Let 
H be a graph on the vertex set {1, . . . ,h}, and let G be a graph on the vertex set 
{1, . . . , n}, where n > h. Let W C V(G) with \W\ = h, and let r\y denote the 
least element in W. Following |T2j, we say that H appears at W in G if (a) the 
increasing bijection from {1, . . . , h} to W gives an isomorphism between H and the 
induced subgraph G[W] of G; and (b) there is exactly one edge in G between W 
and the rest of G, and this edge is incident with the root rw- 

Let an(G) be the number of appearances of H in G, that is, the number of sets 
W C V{G) such that H appears at W in G. Let a be {9e 2 {h + 2))- 1 p h /h\. It is 
shown in |12| that if G n is a random planar graph on n vertices then 

PrW(G„) < an} < e~ an , 

for n large enough. The next result describes more precisely the asymptotic behav- 
ior of the number of appearances of H in random planar graphs. 

Theorem 4. Let H be a fixed rooted connected planar graph with h vertices. Let X n 
denote the number of appearances of H in a random planar graph with n vertices. 
Then X n is asymptotically normal and the mean \i n and variance er^ satisfy 

(1.4) fin - —n, a n - pn, 
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where p — 7 1 and 7 is as in Theorem^ Moreover, for every a < p n jn\ and every 
[3 > p n /n\ we have for n large enough 

i,0 



(1.5) Pr{X„ < an} < — — , Yv{X n > f3n} < 

\x{u)p) \x{u)p 

where x(u) is the solution of 

and u is related to Z , where Z is either a or [3, by the equation 

~u—— = Z. 

x(u) 

Another parameter, the number of 2-connected components in a random con- 
nected planar graph, also follows a normal distribution. 

Theorem 5. Let X n denote the number of blocks (2-connected components) in a 
random connected planar graph with n vertices. Then X n is asymptotically normal 
and the mean p n and variance a„ satisfy 

(1.6) fj, n ~ (n, a 2 n ~ (n, 
where ( w 0.0390518027. 

Next we turn to a different parameter, the number of connected components in 
random planar graphs. 

Theorem 6. Let X n denote the number of connected components in a random 
planar graph with n vertices. Then X n — 1 is distributed asymptotically as a Poisson 
law of parameter v, where v « 0.0374393660. 

The above result is an improvement upon what was known so far. It is shown 
in J2| that Yn is stochastically dominated by 1 + Y, where Y is a Poisson law P(l); 
Theorem [S] shows that in fact Y n is asymptotically 1 + P(v). The following direct 
corollary to Theorem [S] is worth mentioning. 

Corollary 1. (i) The probability that a random planar graph is connected is asymp- 
totically equal to e~ u w 0.9632528217. (ii) The expected number of components in 
a random planar graph is asymptotically equal to 1 + v « 1.0374393660. 

Our last result is the following. Let A be a family of connected planar graphs, 
and let A(x) = A n x n /n\ be the corresponding generating function. Assume that 
the radius of convergence of A(x) is strictly larger than p = 7 -1 , the radius of 
convergence of C(x); this is equivalent to saying that A is exponentially smaller 
than the family C of all connected planar graphs. 

Theorem 7. Assume A is a family of connected planar graphs that satisfies the 
previous condition, and let X n denote the number of connected components that 
belong to A in a random planar graph with n vertices. Then X n is distributed 
asymptotically as a Poisson law of parameter A(p) . 

If we take A as the family of graphs isomorphic to a fixed connected planar 
graph H with n vertices, then 

n\ x n x n 
W = |Aut(#)| ' 7*T = |Aut(JT)|' 
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where Aut(iJ) is the group of automorphisms of H. In particular, if H is a single 
vertex, we obtain that the number of isolated vertices in a random planar graph 
tends to a Poisson law P(p) = P(7 _1 ). This proves a conjecture by McDiarmid, 
Steger and Welsh J2] . 

As a different application of Theorem we have the following. Recall that B(x) 
is the generating function of 2-connected planar graphs. 

Corollary 2. Let X n denote the number of connected components which are 2- 
connected in a random planar graph with n vertices. Then X n tend to a Poisson 
law of parameter B(p) w 0.0006837025. 

We wish to emphasize that the approach that eventually has led to the enumer- 
ation of planar graphs has a long history. Whitney's theorem [21] guarantees that 
a 3-connected graph has a unique embedding in the sphere; hence the problem of 
counting 3-connected graphs is in essence equivalent to counting 3-connected maps 
(planar graphs with a specific embedding) . This last problem was solved by Mullin 
and Schellenberg |13j using the approach developed by Tutte in his seminal papers 
on counting maps (see, for instance, JH|)- The next piece is due to Tutte \nj\: a 
2-connected graph decomposes uniquely into 3-connected "components" . Tutte's 
decomposition implies equations connecting the generating functions of 3-connected 
and 2-connected planar graphs, which were obtained by Walsh [20], using the re- 
sults of Trakhtcnbrot ^Hj- This in turn was used by Bender, Gao and Wormald 
to solve the problem of counting 2-connected planar graphs; their work is most 
relevant to us and is in fact the starting point of our research. Finally, the decom- 
position of connected graphs into 2-connected components, and the decomposition 
of arbitrary graphs into connected components, imply equations connecting the cor- 
responding generating functions. Analytic methods, together with a certain amount 
of algebraic manipulation, become then the main ingredients in our solution. 

Acknowledgements. We are grateful to Philippe Flajolet for his encouragement 
and useful discussions during our research; to Eric Fusy for his help in deriving 
large deviation estimates and in simplifying the final expressions in Lemma and 
to Dominic Welsh for giving us access to an early version of ^2] ■ Discussions with 
Manuel Bodirsky and Mihyun Kang are also acknowledged. 

2. Preliminaries 

In this section and in the rest of the paper we use the language and basic results 
of Analytic Combinatorics, as in the forthcoming book of Flajolet and Sedgewick 0. 
For the sake of completeness, we state the main results we use in this paper (Corol- 
lary VL1, Theorems IX.10 and IX. 13 in gj). 

Proposition 1 (Transfer Theorem; simplified version). Assume that f{z) is ana- 
lytic in a domain A = A(cj>, R), where R > 1, < 4> < 7r /2 and 

A(<p,R) = {z : z±l,\z\<R, \Arg(z-l)\>4>}. 

If as z — > 1 in A, 

f(z)~(l-z)- a 

then 
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Proposition 2 (Quasi-Powers Theorem; algebraic singularities). Let f(z,u) be a 
bivariate function that is bivariate analytic at (0, 0) with nonnegative coefficients 
there. Assume that f admits inT> = {\z\ < r} x {\u — 1| < e}, for some r > and 
e > 0, the representation 

f(z, u) = A(z, u) + B(z, u)C(z, u)~ a , 

where A, B and C are analytic in T> such that C'(z, 1) — has a unique simple root 
p < r in \z\ < r and B(p, 1) =/= 0. Moreover, neither d z C(p, 1) nor d u C(p, 1) are 0, 
so there exists a nonconstant p(u) analytic at u = 1 such that C(p(u),u) = and 
p = p(l). Finally, p{u) is such that 

P "(i) P '(i) , f P '(iy " 



{2 ' 1] p(i) p(i) ' Vp(i) 

is non-zero. Then the random variable with probability generating function 

,_ [z n ]f(z t u) 

M/(z,i) 

converges in distribution to a Gaussian variable. The mean p n and the standard 
deviation o~ n converge asymptotically to pn and a^/n, where p is —p'(l)/p(l) and 
u is given by $2.1]) . 

Proposition 3 (Local Limit Law; simplified version). Let f(x,u) satisfy the con- 
ditions of Proposition^ Lf p(u) attains uniquely its minimum on the circle \u\ = 1 
at u = 1, then the sequence p n i is asymptotically (\f2itn ct) _1 . 

Now we discuss the generating functions that appear in this paper. Recall that 
<?„, c„ and b n denote, respectively, the number of planar graphs, connected planar 
graphs, and 2-connected planar graphs on n vertices. The corresponding exponen- 
tial generating functions are related as follows. 

Lemma 1. The series G{x), C{x) and B(x) satisfy the following equations: 

G(x) = exp(C(x)), xC'(x) = xexp (B'(xC'(x))) , 
where C'{x) and B'(x) are derivatives with respect to x. 

Proof. The first equation is standard, given the fact that a planar graph is a set of 
connected planar graphs, and the set construction in labelled structures corresponds 
to taking the exponential of the corresponding exponential generating function. 

The second equation follows from a standard argument on the decomposition of 
a connected graph into 2-connected components. Take a connected graph rooted 
at a vertex w; hence the generating function xC'(x). Now v belongs to a set of 
2-connected components (including single edges), each of them rooted at vertex v; 
hence the term exp(£>'). Finally, in each of the 2-connected components, replace 
every vertex by a rooted connected graph; this explains the substitution B' \xC' \x)). 
Details can be found, for instance, in [111 p. 10]. □ 

Let b n ^ q be the number of 2-connected planar graphs with n vertices and q edges, 
and let 

Ex n 
bn,qV q —r 

be the corresponding bivariate generating function. Notice that B(x, 1) = B(x). 
The generating functions C(x,y) and G(x,y) are defined analogously. Since the 
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parameter "number of edges" is additive under taking connected and 2-connected 
components, the previous lemma can be extended as follows. 

Lemma 2. The series G(x,y), C(x,y) and B(x,y) satisfy the following equations: 

d ( d d 

G{x,y) = exp(C(x,y)), x—C(x,y) = xexp l—B(x—C(x,y),y) 

In the remaining of the section we recall the necessary results from pQ. Define 
the series M(x, y) by means of the expression 

M m{x , s) = xV (-±- +t l- -i- (1 .r + ^ ( ';r )2 

\1 + xy l + y (1 + U + V) A 

where U(x,y) and V(x,y) are algebraic functions given by 
(2.3) U = xy(l + V)\ V = y(l + U) 2 . 

In the next result and in the rest of the paper, all logarithms are natural. 
Lemma 3 (Bender et al. PQ). We have 

dB(x,y) _ x 2 f l + D(x,y Y 
[ ' ' dy 2 V 1 + V 

where D = D(x,y) is defined implicitly by D(x,0) = and 

(2.5) M^-Ul±R) + JE*-= . 
v 1 2x 2 D 8 V 1 + y J 1 + xD 

Moreover, the coefficients of D(x,y) are nonnegative. 



There is a small modification in equation (|2.4|l with respect to pQ. We must 
consider the graph consisting of a single edge as being 2-connected, otherwise Lem- 
mas ^ an d [3 would not hold. Hence the term of lowest degree in the series B(x,y) 
is yx 2 /2. 

Let us comment on the previous equations. The algebraic generating func- 
tion M corresponds to (rooted) 3-connected planar maps. The decomposition of a 
2-connected graph into 3-connected components implies equations 1)2. 4|l and l|2.5|l . 
The generating function D(x, y) is that of planar networks, which are special graphs 
with two distinguished vertices. 

We define the following functions of the complex variable t. The appendix con- 
tains additional functions that arc introduced later. 



(1 + 3t)(l - t) 5 

m 3 



(1 + 2t) f t 2 {l~t)(\8 + 36t + 5t 2 ) \ 

(1 + 3t)(l -t) CXP V 2(3 + t)(l + 2t)(l + 3t)V 
a = 144 + 592i + 664i 2 + 135t 3 + 6t 4 - 5i 5 
(3 = 3t(l + f)(400 + 1808i + 2527i 2 + 1155i 3 + 237i 4 + m 5 ) 
3t 2 



Do 



D 



2 



(l-t)(l + 3t) 
48i 2 (l + t)(l + 2f) 2 (18 + 6t + t 2 ) 



(1 + 3i)/3 

D 3 = 384:t 3 {l + t) 2 (l + 2t) 2 (3 + t) 2 a 3/2 p- 5/2 
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Let us notice a slight change in terminology: functions £ and Y are denoted, re- 
spectively, £0 and yo in Dj : also, we correct a typo, namely a t factor that was 
missing in the expression for D 2 . 

A key fact is that for y in a suitable small neighborhood of 1, the equation 
Y(t) = y has a unique solution in t = t(y). Then define 

(2-6) R(y) =£(%)). 

In the next lemma, Di(y) stands for Di(t(y)). This applies too to functions Bi(y) 
and Ci(y) that we introduce later in the paper. 

Lemma 4 (Bender et al. p^). For fixed y in a small neighborhood of 1, R(y) is 
the unique dominant singularity of D(x,y). Moreover, D(x,y) has a branch-point 
at R{y), and the singular expansion at R{y) is of the form 

D(x, y) = D (y) + D 2 (y)X 2 + D 3 (y)X 3 + 0(X 4 ), 

where X = \/l — x/R{y) and the Di{y) are as before. 

The previous lemma is the key result used in pQ to prove the estimate 

b n ~b-n- 7 / 2 R- n nl, 

where 6 is a constant and R = i?(l) w 0.0381910976. 



3. Analysis of B(x,y) 
From equation (|2.4fl . it follows that 
, , , , x 2 f v 1 + D(x.t) , 

(3.1) b ^ = tJ -^tt 2 ^ 

Our goal is to obtain an expression for B(x,y) as a function of x,y and D(x,y) 
that, although more complex, does not contain an integral. Recall that the algebraic 
function U is defined in (|2.3|l . and D is defined in Lemma [21 

Lemma 5. Let W(x,z) = z(l + U(x 1 z)). The generating function B(x,y) of 2- 
connected planar graphs admits the following expression as a formal power series: 



(3.2) 

where 

and 

f3 2 (x,z,w) 



B(x, y)=f3 (x, y, D(x, y), W (x, D(x, y))) 



X X 

(3(x,y,z,w) = —j3x{x,y,z) - -f3 2 (x,z,w), 



z{6x-2 + xz) f 1 + y\ 1 °g( 1 + ^) 1 log(l + a;z) 

4x +[ +Z > g \l + z) 2 + 2x 2 : 

2(1 + x)(l + w)(z + w 2 ) + 3(to - z) 1 



2(1 + wf 



I -Ax 
2x 



log(l + w) 



1-Ax+ 2x 2 
Ax 



log 



— log(l + xz + xw + xw 
2x 

1 — x + xz H — xw + xw 2 
(1 - x)(z + w 2 + 1 + w) 
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Proof. From equation 1|3.1|1 we obtain 
We integrate by parts and obtain 

jf TTT d< = log(1 + y)D{x > y) ~ Jo log(1 + Q^at^ dt 

From now on a; is a fixed value. Now notice that from 1)2. 5JI it follows that 

M(x,u) xu 2 
2x 2 u 1 + xu 



(u) = -1 + (1 + u)exp - 



is an inverse of D(x,y), in the sense that 4>(D(x,y)) = y. In the last integral we 
change variables s = D(x, t), so that t — 4>(s). Then 

dD(x, i) f D(x ' y) 



1 + xs 

D{x ' v) M{x, ^ 



c/.s 



2x 2 i 



ds. 



The first integral has a simple primitive and we are left with an integral involving 
M(x, y). Summing up we have 

(3.3) B{x,y)=Q{x,y,D{x,y)) + - ^^-ds, 

4 7o s 

where is the elementary function 



0(1, y, z) = — + -z 2 + (l + z) log ^ ~2 Z+ 2 bg(1 + XZ) ' 
Now we concentrate on the last integral. From (|2.2|l and (|2.3(l it follows that 
^ds = xf D ^t U)2 Lds, 



where [/ and V are considered as functions of x and s, and where for simplicity we 
write D = D(x, y) from now on. 

From the definition W(x, s) = s(l + U(x, s)), we obtain that 

{1 + U) 2 U W-s 
(1 + U + Vf ~ W{l + W) 3 ' 
Since W satisfies the equation 

xs 2 + (1 + 2xW 2 )s + W{xW 3 - 1) = 0, 
the functional inverse of W(x, s) with respect to the second variable is equal to 



_ 2 _ 1- Vl + 4xf + 4xt 2 
ld j 2a; 
where we use t to denote the new variable. 
It follows that 

f D W-s rW{x,D) (Q-\-2xt- 2xt 2 ) (2Qt - 2* - 1) 

as = / — - r^-^ at, 



o ^(1 + 1^)3 ' J Q 2xt(l+t) 3 ( 
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where for simplicity we write 



(3.5) Q{x,t) = y/l + Axt + Axt 2 . 

The last integral can be solved explicitly with the help of a computer algebra system 
such as Maple, and we obtain as a primitive the function 

1 - 2{t + Ax + Axt) 1 + 2x(l + t)~ / A 1 + 2(t - x - tx) \ 

-Q 3 + 2 + 4xt + -V- — — - Q + 



4x(l + t) 2 2x(l + t) 2 ^ V 4x(l + t) 2 

2x 2 - Ax + 1 {Q + (l-2x- 2xt) \ 1 , , , 1 - 4x 

Finally we have to replace t for W(x, D) in the previous equation. The expression 
(|3.4|) and equation (|3.5fl imply that 

Q(x, W(x, D)) = l + 2x(D + W{x, D) 2 ). 

Hence when replacing t for W(x, D) we obtain an expression in x, D and W(x, D) 
that is free of square roots. A routine computation, combined with the intermediate 
equation lj3.3|l . gives the final expression for B(x,y) as claimed. □ 

The function (3 in the previous lemma looks like the primitive of a rational 
function. This can be explained as follows (we are grateful to P. Flajolet for this 
observation). The algebraic equation satisfied by U (here x is considered as a 
parameter) is 

u-xy(l + + uff = 0. 
It can be checked (for instance, using the Maple package algcurves), that this 
equation in u and y defines a rational curve, that is a curve of genus zero, and so it 
admits a rational parametrization (u(t),y(t)). Now an integral J R(s,U(x, s)) ds, 
where R is a rational function, becomes the integral of a rational function after 
the change of variables s — y(t). In particular, this applies to the integral in 
equation l|3.3(l . 

The former lemma can be used to obtain the singular expansion of B(x, y). The 
function R(y) is defined in (|2.6|l and Bq, B-2, B4, B$ are analytic functions of y given 
in the appendix. Again Bi(y) stands for Bi(t), where t is the unique solution of 
Y(t) = y in a neighborhood of 1. 

Lemma 6. For fixed y in a small neighborhood of 1, the dominant singularity of 
B(x, y) is equal to R(y). The singular expansion at R(y) is of the form 

B(x, y) = B (y) + B 2 (y)X 2 + B 4 (y)X 4 + B 5 (y)X 5 + 0(X 6 ), 



where X = yl — x/R(y), and the B{ are analytic functions in a neighborhood of 1. 

Proof. Consider the expression for B(x,y) in Lemma [5] as a function of x,y and 
D(x,y). A simple analysis shows that, for y close to 1, the only singularities come 
from the singularities of D(x,y), hence the first claim follows. 

For the second assertion, substitute the singular expansions of D(x, y) and 
U(x, D(x.yj) (taken, respectively, from Lemma 0| and the appendix) for D(x,y) 
and U(x,D(x,y)) in (O (recall that W = z{l + U)). Next set x = £(t)(l - X 2 ) 
and y = Y(t) as functions of t, and expand the resulting expression. That the 
coefficients Bi are as claimed in the appendix is a tedious but routine computation 
that we performed with the help of Maple. In particular, the coefficients of X and 
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X 3 vanish identically in y (or in t). The Bi are analytic since they are elementary 
functions of the Di . □ 

4. Asymptotic estimates 

In order to prove Theorem ^ first we need to locate the dominant singularity 
p = 7 _1 of G(x). Since G(x) — exp(C(a:)), the functions G{x) and C{x) have the 
same singularities; hence from now on we concentrate on C{x). 

We rewrite the second equation in Lemma ^ as 

(4.1) F{x) = xex-p(B'(F(x))), 

where F(x) = xC'(x). Notice that the singularities of B'(x) and F(x) are the same, 
respectively, as those of B(x) and C(x). From j|4.1(l it follows that 

(4.2) ff>(u) = ue- B ' {u] 

is the functional inverse of F(x). The dominant singularity of if) is the same as 
that of B(x), which according to Lemma El is equal to R = B(l). In order to 
determine the dominant singularity p of F(x), we have to decide which of the 
following possibilities hold; see Proposition IV. 4 in |H] for an explanation. 

(1) There exists r € (0,B) (necessarily unique) such that tp'( T ) = 0- Then ip 
ceases to be invertible at r and p — ip{r). 

(2) We have f(u) ^ for all u S (0, R). Then p = ip(R). 

The condition iP'(t) = is equivalent to B"(t) = 1/t. Since B"(u) is increasing 
(the series B(u) has positive coefficients) and 1/u is decreasing, we are in case (2) 
if and only if B"(R) < 1/R. Next we show that this is the case. 

Claim 1. Let R be as before the radius of convergence of B(x). Then B"(R) < l/R. 

Proof. Lemma implies that B"(R) = 2£? 4 /i? 2 (see l|4.3|l below). Hence the in- 
equality becomes 2B 4 < R. It holds because R w 0.0381 and B 4 w 0.000767. □ 

Let us remark that in a related problem, counting series-parallel graphs, a very 
similar situation appears but the analogous ip function does have a maximum in its 
domain of definition |2). 

We are now ready for the main result. 

Proof of Theorem^ As we have seen in the previous claim, the dominant singularity 
of F(x) is at p = ip(R). In order to obtain the singular expansion of F{x) at p, we 
have to invert the singular expansion of ip(u) at R. 

The expansion of B'(x) follows directly by differentiating the one in LemmaEl 

(4.3) B'(x) = -i (b 2 + 2B 4 X 2 + \B 5 X^j + G(X 4 ). 
Because of ip( x ) = xexp(—B'(x)), by functional composition we obtain 

i>{x) = Re B */ R (l + l) X- + 5 -^X^ + 0{X% 

Since we are inverting at the singularity, F(x) also has a singular expansion of 
square-root type 

F{x) =F + F X X + F 2 X 2 + F 3 X 3 + 0(X 4 ), 
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with the difference that now X = y/l — x/p. Given that F(x) and ip(x) are func- 
tional inverses, the Fi can be found by indeterminate coefficients, and they turn 
out to be, in terms of R and the Bi, 

(4.4) F = R, Fi = 0, F 2 = 2B ^_ R , F 3 = -|b 5 (1 - 2B 4 /R)- 5 / 2 . 

The singular expansion of C{x) is obtained by integrating C'(x) = F(x)/x, and 
one gets 

(4.5) C{x) = C + C 2 X 2 + C A X A + C 5 X 5 + 0(X 6 ). 

The C, except Co, are computed easily in terms of the Fi in equation 14. 411 . and 
they turn out to be 

Fn + F 2 2 

(4.6) C 2 = -F , C 4 = -±-°±^, C 5 =--F 3 . 

2 5 

By singularity analysis (Proposition 0, we obtain the estimate 

c„ ~ c • n~ 7/2 p~ n n\, 

where c = C 5 /r(-5/2). 

However, the coefficient Co = C(p) is indeterminate after the integration of 
F(x)/x, and is needed later. To compute it, we start by integrating by parts 

C(x)= [ ^-ds = F(x)logx- [ F'(s)logsds. 
Jo s Jo 

We change variables t = F(s), so that s = ip(t) = te~ B W, and the last integral 
becomes 

logip(t)dt = / (log* - B' '(*)) dt = F(x) log F(x) - F(x) - B(F(x)). 
'0 Jo 

Hence 

C(x) = log x - F(x) logF(x) + if (a:) + B(F(x)). 
Taking into account that F(p) = R and B(R) = Bo, we get 

Co =C(p) = R log p - R log R + R + Bo. 
A simple computation shows that, equivalently, 

(4.7) C = R + B + B 2 . 

The final step is simpler since G(x) — e ctyX \ We apply the exponential function 
to 14. 5f) and obtain the singular expansion 

(4.8) G{x) = e c ° (l + C 2 X 2 + (C 4 + ^C 2 2 )A 4 + C 5 X 5 ^j + 0(X 6 ), 

where again X = yl — x/p. Again by singularity analysis, we obtain the estimate 

g n ~ 9 ■ n- 7/2 p- n nl, 

where g = e c °c. Finally, since p = ip(R) = Rer B ' (R) and B'(R) = -B 2 /R, we get 

p = Re B ^ R , 1 = P - X = -e- B -' R . 

R 

The constants c, g and p can be found using the known value of R and the expres- 
sions for the Bi in the appendix; the approximate values in the statement have 
been computed using these expressions. 
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Notice that the probability that a random planar graph is connected is equal to 
Cn/g-n ~ c/g — e~ c °. This result reappears later in Theorem^! □ 

5. Gaussian limit laws 

The proofs in this section are based on bivariate singular expansions and per- 
turbation of singularities. To simplify the notation, in this section we denote by 
f'{x,y) the derivative of a bivariate function with respect to x. 

Proof of Theorem^ We rewrite the second equation in Lemma |21 as 

(5.1) F(x,y) = xexp(B'(F(x,y),y)), 
where F(x,y) = xC'(x,y). It follows that, for y fixed, 

(5.2) ip{u,y) = ue~ B ' {u ' v) 

is the functional inverse of F(x,y). 

We know from the previous section that ip'(u,y) does not vanish for y = 1 and 
u 6 (0, R), and that p — ip(R) is the dominant singularity of F(x). Hence by 
continuity the same is true for y close to 1, and the dominant singularity of F(x, y) 
is at 

(5.3) P (y) = 4>{R{y), y) = R{y)e- B '^^\ 

Given the analytic expressions for the functions involved, the univariate singular 
expansion of i/j(x) extends to an expansion of i/)(x, y) for y fixed. The same is true 
then for F(x,y) and C(x,y), and we obtain a bivariate expansion 

C(x, y) = C (y) + C 2 {y)X 2 + C i {y)X i + C 5 (y)X 5 + 0(X 6 ), 

where the Ci(y) are analytic functions, and now X = v/l — x/p(y). 

Then Proposition|2]implies a limit normal law for the number of edges in random 
connected planar graphs, with expectation and variance linear in n. The constants 
K and A in the statement of Theorem |21 are given by 

\ p"(l) p'(l) fp'(l)\ 2 
K P(l)' P(l) \P(1)J ' 

where p'(y) — dp(y)/dy. Since G(x,y) and C(x,y) have the same dominant sin- 
gularities p(y), the previous statement also holds for arbitrary planar graphs, with 
the same values of k and A. 

In order to determine the parameters exactly, we need only an explicit expression 
for p(y). The expansion l|4.3ll extends to an expansion of B'(x,y), whose constant 
term is B'(R(y),y) — —B2(y)/R(y). Hence from l|5.3|) it follows that 

p(y) = R(y)vq>(B 2 (y)/R(y)). 

The appendix contains an explicit expression for p{y) = q(t) as a function of t. 
The necessary derivatives are computed as p'{y) — q'(t)/Y'(t), and the same goes 
for p"{y). The approximate values in the statement have been computed in this 
way. □ 
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Proof of Theorem\^ Consider the generating function C u (x,y) = C(x,uy), where 
u is a fixed constant. In this situation the singularity p u {y) of C u is given by p{uy), 
and the associated probabilities k of C u are 

u _ [y k ][x n ]C u (x,y) _ u k c n , k 
[ ' ' Pn > k [x n ]C a (x,l) n\[x n ]C(x,u)' 

In order to apply Proposition [3] to C u we need to know the singularities of C(x,y) 
when y is away from 1. The following claim extends Claim ^ an d shows that the 
bivariate singularity expansions given in the proof of Theorem |3 hold for every y. 

Claim 2. Let R(y) be the radius of convergence of B(x,y) for y fixed. Then 
B"(R(y),y)<l/R(y). 

Proof. As in the proof of Claim it is enough to show that 2B 4 {y) < R(y) for 
y G (0,oo); equivalently, that 2i3 4 (i) < £(t) for t £ (0, 1). We bound the logarithm 
that appears in the expression for B4 (see the appendix) as 



VT+2I/ ~ VT+21 

Let B4 be the function obtained by substituting the logarithm in B4 for the right- 
hand side in the previous inequality. Then it is enough to show that 

2£ 4 (*)<£(i) forte (0,1). 

Since both B4 and £ are rational functions, the problem reduces to showing that 
a certain polynomial (in fact, of degree 20) with integer coefficients has no root 
in (0, 1). We have checked that this is indeed the case using Maple. □ 

Another requirement is that p(z) attains uniquely its minimum on \z\ = u at 
z = u. Suppose it exists w 7^ u with \w\ = u such that |p(u>)| < p(u). It follows 
from (|5.3I) that R{z) is equal to F(p(z), z), and since F(x,y) has non-negative 
coefficients, |-R(u>)l = \F(p(w),w)\ < F(p(u),u) — R(u). However, this contradicts 
the fact that R(z) attains uniquely its minimum on \z\ = u at z = u, as shown 
in ^ Lemma 3]. 

Now Proposition |21 applied to C u yields 



1 



(5.5) 

where p{u) are o~(u) are given by 

/ x P'ui 1 ) U P'( U ) 

/Ml) p{u) 

f s2_ PuM gUl) , (PuWV _ 2 P» p'(u) ( p'{u) 

^ Pu(l) Pu(l) + \pu(l)J - U P(U) U p(u) + { U p(u). 

Theorcm[3]follows by combining equations (|5.4|l and (|5.5() for k = [p(u)n\ and using 
the asymptotic expression of [x n ]C(x,y) for y = u. The value p is constrained to 
the interval (1, 3) since lim u ^o p(u) — 1 and liniu^oo p(u) = 3. 

Proof of Theorem^ Let us recall Equation 1|4.1|) 

F(x) = xexp(B'(F(x))), 
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where F(x) = xC'{x) is the generating function of rooted connected planar graphs. 
In order to mark appearances of H, we have to look at the root r of a rooted 
connected graph G, and the blocks to which it belongs; recall this is encoded in the 
term exp(B' (F (x))) . We are interested in the blocks which are equal to a single 
edge rv, and within these blocks to the situation where vertex v is substituted 
by a copy oi H. In this case we mark an appearance of H with the secondary 
variable y. If we let f(x,y) be the corresponding generating function, then the 
previous discussion translates into the equation 

(5.6) f(x, y) = x exp ^B'(f(x, y)) + (y - 1)^ 

Notice that from the definition of appearances we do not need to take into account 
the automorphisms of H; if a copy of H is substituted for vertex v, there is only 
one way to do it once the labels are selected, hence the term x h /h\. 

In fact, f(x, y) is not the exact counting series, since it does not take into account 
the possibility that the root r belongs to a copy of H that appears in G. This can 
be accounted for as follows. The generating function of rooted connected graphs 
where the root belongs to an appearance of H is f(x, y)x h /(h — 1)!, since the root 
can appear in any of the h vertices of H . Hence the generating function that counts 
exactly all appearances is 

x h 

9(x, y) = f(x, y) + {y- 1 ) ( h _ 1 y/ ( x > v)- 

Since f(x, y) and g(x, y) have the same dominant singularity for any fixed y it does 
not matter which one we choose for singularity analysis; hence in the rest of the 
proof we work with f(x,y), defined through (|5.6() . 
Equation (|5.6|) can be rewritten as 

f{x, y) = C(x, y) exp (B'(f(x, y))) , 

where ((x,y) — xexp((y — l)x h /h\). Comparing the previous equation with 14. If) , 
it follows that 

f(x,y) = F(((x,y)). 

Given that p is the dominant singularity of F(x), the dominant singularity of f(x, y) 
for fixed y is the smallest value r(y) satisfying 

(5-7) C(r(j/), V) = r(y) exp ((y - l) 1 ^) = P- 

Clearly r(l) = p. In order to compute r'(y), we differentiate l|5.7[l . set y = 1, and 
obtain r'(l) = -p h+1 /h\. To compute r"(l) we differentiate again and, after a 
simple computation, we get 

r'(l) p h t"(1) r'(l) /r'(l) 



r(l) h\ ' t(1) r(l) ' \t(1) J P ' 

From the singular expansion of F(x) at p, we derive a corresponding bivariate 
singular expansion of f(x,y) at r(y), and again a normal limit law follows from 
Proposition As in the previous proof, a large deviation estimate also follows, 
and from this we obtain the bounds in (jl.5JI : the details are omitted to avoid 
repetition. 
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Proof of Theorem^ The proof is similar to the previous proofs, and so we omit some 
details. The generating function Ci(x,y) of connected planar graphs according to 
the number of vertices and blocks satisfies the equation 

xC[ (x,y)=x exp (y B' (xC[ (x, y))) , 

where B(x) is the univariate generating function of 2-connected planar graphs. 
Let Fi(x,y) — xC[(x,y). Then, for y fixed, 

is the functional inverse of F\(x, y). The dominant singularity of ipi(u,y) is at R, 
which in this case is independent of y, and the dominant singularity of Fi(x,y) is 
at 

Pl (y) = MR,y) = R<>- yB ' {R) - 

Again we have bivariate singular expansions whose coefficients are analytic func- 
tions of y, and the quasi-powers theorem implies asymptotic normality of the param- 
eter. The asymptotic expressions for the expected value and variance are obtained 
as before, but in this case the computations are particularly easy, since 

p[(y) = - P i(y)B'(R). 

We know that p = tp(R) = Re^ B ' {R \ hence 

C = = B'(R) = \og(R/p) « 0.03905180273. 

Pil 1 ) 



A similar computation gives 

= B'{R) = (. □ 



p'l(i) Pi(i) , (pW 2 



pi(i) VMi). 

6. POISSON LIMIT LAWS 

As opposed to the proofs in the previous section, to prove Theorems El and 
univariate asymptotics is enough. 

Proof of Theorem^ Let v = C(p) = Co, the evaluation of C(x) at its dominant 
singularity. For fixed k, the generating function of planar graphs with exactly k 
connected components is 

liO(x)K 

For fixed k we have 

[x k ]C{x) k - kC%- l [x n ]C{x). 

Hence the probability that a random planar graphs has exactly k components is 
asymptotically 

[x n ]C(x)/k\ kCt 1 Co= v k - x v 
\x n \G(x) kl (fc-1)! ' 

as was to be proved. □ 

Proof of Theorem The proof is similar to the previous one. The generating 
function of planar graphs with no component belonging to A is exp(C(a;) — A(x)). 
Hence the generating function of planar graphs with exactly k components in A is 

^A(x) k exp(C(x) - A(x)) = ^A(x) k e - A{x) G(x). 
kl k\ 
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The same kind of simple calculation as before gives that the probability that a 
random planar graphs has exactly k components in A is asymptotically 

k\ 

This finishes the proof of the theorem. □ 



7. Concluding remarks 

We have found a solution to the problem of counting labelled planar graphs; 
however, counting unlabelled planar graphs appears to be much more difficult. If 
u n is the number of unlabelled planar graphs on n vertices, then it is known that 
the following limit exists 

7 U = lira {u n ) 1/n , 

n — >oc 

and that 7 < 7„, where 7 is as in Theorem ^ (see ^U). The reason for the strict 
inequality 7 < 7„ is that, contrary to what happens for unrestricted graphs, a planar 
graph has with high probability an exponential number of automorphisms |12) . 

The best upper bound obtained so far is 7^ < 30.06. This is proved in [3] by 
showing that an unlabelled planar graph with n vertices can be encoded with 4.91n 
bits. On the other hand, our determination of 7 provides a lower bound on j u , and 
shows that at least 4.76 w log 2 7 bits per vertex are needed. 

We believe that to determine j u exactly is a very hard problem, not to speak 
of determining the subexponential behavior of u n . The reason is that the equa- 
tions connecting the generating functions of labelled planar graphs with different 
connectivity requirements, do not hold anymore in the unlabelled case. 

A related problem is to estimate the number of planar graphs with a given 
number of edges. In |2j it is proved that an unlabelled planar graph with m edges 
can be encoded with 2.82m bits. We can show that at least 2.59 bits per edge are 
needed, as follows. 

The coefficient of y rn in G(l, y) is equal to 

1 V - ^ 9n,m 

TV. 

n 

where g n ,m is the number of labelled planar graphs with n vertices and m edges. 
Since a graph on n vertices has at most nl automorphisms, the number of unlabelled 
planar graphs with m edges is at least h m . 

The exponential growth of the h m is determined by the smallest singularity r 
of G(l,y). Since the smallest singularity of G(x,y) for fixed y is p(y), as given 
in H5.3|l . it follows that r is the smallest solution to p(r) = 1. It can be computed 
exactly with the expressions in the appendix and it turns out that 

lim {h m ) 1/m = t m 6.03 « 2 2 59 . 

n — >oo 

Finally, let us mention that the explicit expressions we have obtained for the 
generating functions of labelled planar graphs have been applied to the design of 
very efficient algorithms for generating random planar graphs uniformly . 
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Appendix 

Here we list the functions So, B2, -B4, -B5 that have been used in the previous 
sections, as functions of t. As has been explained already, to become functions 
of y near 1, they must be evaluated at the unique solution of Y(t) = y. For 
completeness, we also write down the function p. 



Bn = 



(3f- l) 2 (l + t) 6 log(l + t) _ (3t 4 -16* 3 + 6* 2 -l)log(l + 3£) 
512 1 6 32t3 
(l + 3t) 2 (l-t) 6 log(l + 2t) 1 . , 1, . , 3, ,,„, 

(217 1 6 + 920 1 5 + 972 1 4 + 1436 t 3 + 205 1 2 - 172 1 + 6) (1 - tf 
2048i 4 (l + 3i) (3 + i) 



B 2 



(1 - t) 3 (3 1 - 1) (1 + 3 1) (1 + tf log (1 + t) 
256 i 6 

(l-t) 3 (l + 3t)log(l + 3t) (l + 3t) 2 (l-t) 6 log(l + 2t) 



+ 



32 i 3 512 
(1 - tf (185 i 4 + 698 t 3 - 217 i 2 - 160 1 + 6) 
1024t 4 (l + 3t) (3 + i) 



= TT^-TZ r 



512 i 6 2048t 4 (3 + t)Q 



V3 (1-t) 6 / S \ 5/2 



Ss ~" 90 (i+t) 3 / 2 v^y 



where 



P = -2400 + 57952 1 + 303862 i 2 + 466546 t 3 + 264775 i 4 + 76679 1 5 + 1 1495 t 6 + 739 t 7 
Q = 400 + 1808t+ 2527 i 2 + 1155t 3 + 237 t 4 + 17t 5 
S = 144 + 592t + 664i 2 + 135t 3 + 6< 4 -5t 5 

p = ^-VT+3t(l-tft- 3 cxp(A) 7 
16 

where 

A log(l+t)(3t-l)(l+t) 3 | log(l + 2t)(l + 3t)(l-Q 3 



+ 



16 1 3 32 1 3 

(1 - t)(185t 4 + 698i 3 - 217t 2 - 160* + 6) 
64t(l + 3t) 2 (3 + t) 
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The approximate values of the univariate constants Bi = Bi(t), where t — 0.6263716633 
is the unique solution of Y(t) = 1, are 

B a = 0.7396995711 • 10~ 3 , B 2 = -0.1491431215 • 10~ 2 , 
B A = 0.7671782851 • 10~ 3 , B 5 = -0.3501857790 • 10~ 5 . 

Finally we include the singular expansion of U(x, D(x,y)) at the dominant sin- 
gularity R(y): 

U(x, D(x, y)) = U (y) + U x {y)X + U 2 {y)X 2 + 0(X 3 ), 



where X = \J\ — x/R(y) and the Ui are, as functions of t, given by 
1 

3? 



Uq = 
Ui = 



4(1 + 3t) 2 (— 5t 5 + 6i 4 + 135i 3 + 664i 2 + 592i + 144) N 1/2 



27i 3 (l + t)Q 



U 2 = 



2(l + 3t)T 



27t 2 {l + t) 2 Q 2 
where Q is as before and 

T = 691t 12 + 10112i n + 98693< 10 + 719346i 9 + 3723625i 8 + 13180580t 7 + 31133003t 6 
+47691938i 5 + 47354348t 4 + 30156200i 3 + 11835336t 2 + 2596736£ + 243072 
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